COMPUTING MODULAR EQUATIONS FOR SHIMURA CURVES 



YIFAN YANG 



ABSTRACT. In the classical setting, the modular equation of level N for the modular 
curve Xo(l) is the polynomial relation satisfied by j(r) and j(Nr), where j(r) is the 
standard elliptic j'-function. In this paper, we will describe a method to compute modular 
equations in the setting of Shimura curves. The main ingredient is the explicit method for 
computing Hecke operators on the spaces of modular forms on Shimura curves developed 
in IB). 



1. Introduction 
Let j(r) be the elliptic j -function. For a positive integer N, consider the set 



T N =SL(2,Z) ( I °)SL(2,Z) 



a b 

c d 



N 

a, b, c, d G Z, ad — be = N, gcd(a, b, c, d) — 1 



The group SL(2, Z) acts on Tjv by multiplication on the left and it can be easily checked 
that the cardinality of SL(2, Z)\IV is finite. Moreover, the multiplication of any element 
in SL(2, Z) on the right of the cosets in SL(2, Z)\rjy simply permutes the cosets. Thus, 
any symmetric sum of j'(7 T )> wnere 7 runs over a complete set of representatives for 
SL(2, Z)\rjv, will be a modular function on SL(2, Z) and therefore can be written as a 
rational function of j(r). In fact, because the only possible pole occurs at the cusp, this 
rational function of j(r) is actually a polynomial. Then the modular polynomial ^^(x,y) 
of level N is defined to be the polynomial in C[a;, y] such that 

7GSL(2,Z)\r N 

In fact, it can be easily proved that this polynomial $>n{x, y) is in Z[x, y\. Since Yq(1) is 
the moduli space of isomorphism classes of elliptic curves over C, the equation $n (x, y) — 
is the modular equation of level N for moduli of elliptic curves over C. That is, if Ei 
and E2 are two elliptic curves over C admitting a cyclic 7V-isogeny between them, then 
their j-invariantsj(Si) and j(E 2 ) satisfy $ n (j (Ei) , j (E2)) — 0. 

Observe that if the Fourier expansion of j{r) is g _1 + 744 + 196884q+ • • • , q = e 2, ™ T , 
then the Fourier expansion of j(Nr) is simply q~ N + 744 + 196884^ + • • • . Thus, in 
principle, to determine §>n(x, y), one just has to compute enough Fourier coefficients for 
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j(r) and solve a system of linear equations. Of course, the main difficulty in practice is 
that the coefficients are gigantic. On the other hand, because the coefficients of <&n{x, y) 
are all integers, one can compute the reduction of $n (x, y) modulo p for a suitable number 
of primes p and then use the Chinese remainder theorem to recover the coefficients. See 
[ 3 1 for the current state of the art in the computation of $n(x, y)- 

In this paper, we shall consider modular equations in the settings of Shimura curves. Let 
B be an indefinite quaternion algebra of discriminant D > 1 over Q and O be a maximal 
order in B. Choose an embedding i : B M(2, R) and let 

r(O) := {t(a) : a e O, N{a) = 1} 

be the image of the norm-one group of O under l, where N(a) denotes the reduced norm 
of a G B. Then the Shimura curve X®(1) is defined to be the quotient space r(0)\H. 
As shown in [9], the Shimura curve X D (1) is the coarse moduli space for isomorphism 
classes of abelian surfaces with quaternionic multiplication (QM) by O. Let Wd denote 
the group of Atkin-Lehner involutions on Xq(1). For our purpose, we will also consider 
the quotient curves of Xq (1) by subgroups W of Wd- 

Now assume that X D (l) has genus and choose a Hauptmodul t(r) for X®(1) so that 
t{r) generates the function field on X®(1). For a positive integer N relatively prime to 
D, pick an element a of norm N inO such that O n (a _1 Oa) is an Eichler order of level 
N. Then the modular polynomial of level N for X D (1) is defined to be the polynomial 
®n( x 'V) °f minimal degree, up to scalars, such that $^(t(r), t(i(a)r)) = 0, which is 
essentially the rational function ^^(x, y) such that 

*£(s,i(T))= II fr-Krr)). 

7 6r(0)\r(o)i(a)r(o) 

Here unlike the case of the modular curve Xo(l), a symmetric sum of t(jr) as 7 runs 
through representatives of T(0)\T(0)i(a)T(0) is not equal to a polynomial of t{r) in 
general. Since X D {1) is the moduli space of abelian surfaces over C with quaternionic 
multiplication, the modular equation of level N relates the moduli of two abelian surfaces 
with QM that have a certain type of isogenies between them. (The precise description 
of the isogeny is a little complicated to be given here. See |7, Appendix A] for details.) 
Modular equations for Atkin-Lehner quotients X® (1)/W are similarly defined. 

When D > 1, the problem of explicitly determining modular equations for Shimura 
curves is significantly more complicated than its classical counterpart. The reasons are 
that Shimura curves do not have cusps and it is difficult to find the Taylor expansion of an 
automorphic function with respect to a local parameter at any given point on the Shimura 
curve. When D and N are both small such that Xq(N) is also of genus 0, it is possible 
to work out an explicit cover X ( f(N) — > X D (1) from the ramification data alone. Then 
from the explicit cover, one can compute the modular equation of level N. This has been 
done in [4] for a limited number of cases. 

Another possible method to compute modular equations for Shimura curves uses the 
Schwarzian differential equation associated to a Hauptmodul i(r). (See Section [2] for a 
review on the notion and properties of Schwarzian differential equation.) The idea is that 
with a properly chosen pair of solutions F\ (t) and F2 (t) of the Schwarzian differential 
equation and a correct positive integer e, the expression (F2(t) / Fi(t)) e can be taken to 
be a local parameter at the point tq of the Shimura curve with t(ro) = 0. Inverting the 
expression, one gets the Taylor expansion of t with respect to the local parameter. If 
somehow one manages to find the Taylor expansion of t(i(a)r) (this can always be done 
when i(t(a)ro) = 0, then by computing enough terms and solving a system of linear 
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equations, one gets the modular equation. However, as far as we know, there does not 
seem to be any paper in the literature that employs this idea to obtain modular equations 
for Shimura curves. (The paper [1 ] did obtain the Taylor expansion for the Hauptmodul in 
the case D = 6. More recently, Voight and Willis [13] developed a method for numerically 
computing Taylor expansions of automorphic forms on Shimura curves.) 

In his Ph.D. thesis |fl2]| . Voight developed a method to compute modular equations for 
Shimura curves associated to quaternion algebras over totally real number fields in the 
cases when the Shimura curves have genus zero and precisely three elliptic points, i.e., 
Shimura curves associated to arithmetic triangle groups. The idea is to use hypergeo- 
metric functions, which are essentially solutions of the Schwarzian differential equations 
mentioned above, to numerically determine coordinates of CM-points and then use the 
existence of canonical models and explicit Shimura reciprocity laws to determine modu- 
lar equations. Even though the equations are obtained numerically, they can be verified 
rigorously by showing that the monodromy group of the branched cover is correct. 

In this paper, we shall present a new method to compute modular equations for Shimura 
curves. Our method also uses Schwarzian differential equations, but relies more heavily on 
the arithmetic side of the theory. Namely, in an earlier work [ 15 ], we showed that the spaces 
of automorphic forms of any given weight on a Shimura curve X®(1) of genus can be 
completely characterized in terms of solutions of the Schwarzian differential equation. We 
lfT31l then devised a method to compute Hecke operators with respect to our basis. The 
Jacquet-Langlands correspondence plays a crucial role in our approach. It turns out that 
our method and results in [ 15] can also be used to compute modular equations for Shimura 
curves. We will describe the procedure in Section[3]and give a detailed example in Section 
E] Using the computer algebra system MAGMA [2|, we have succeeded in determining 
modular equations of prime level up to 19 for Xq (I) /We and those of prime level up to 23 
forX 10 (l)/^ 10 . 

The main difficulty in generalizing our method to general Shimura curves lies at the 
fact that our method requires that a Schwarzian differential equation is known before- 
hand. Because of the problem of the existence of accessory parameters, the determination 
of Schwarzian differential equations usually requires that an explicit cover X^N) — > 
Xq(1) is known, which can be problematic when D is large. (In some sense, what we 
do here is to deduce modular equations of higher levels from that of a given small level.) 
Nonetheless, once an explicit cover of Shimura curves is determined, the combination of 
the methods in [15] and in this paper will yield modular equations for the Shimura curve. 

The rest of the paper is organized as follows. In Section |2] we review the definition 
and properties of Schwarzian differential equations. In Section[3] we describe our method 
to compute modular equations for Shimura curves X® (\) /W , assuming that an explicit 
cover Xq(pq)/W — > Xq{1)/W is known. In Section HI we give a detailed example 
illustrating our method. In Sections and [6] we list part of our computational results for 
the Shimura curves Xq(1)/W 6 and Xq°(1)/Wi . Files in the MAGMA-readable format 
containing all our computational results are available upon request. 



2. Schwarzian differential equations 

In this section, we will review the definition and properties of Schwarzian differential 
equations. In particular, assuming the Shimura curve has genus 0, we will recall the char- 
acterization of the spaces of automorphic forms in terms of solutions of the associated 
Schwarzian differential equation. The method for computing Hecke operators on these 
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spaces is too complicated to describe here. We refer the reader to lfT5l for details. Most 
materials in this section are taken from lfl5l . 

Let X be a Shimura curve. We assume that the associated quaternion algebra is not 
M(2, Q) so that X has no cusps. Let F(t) be a (meromorphic) automorphic form of 
weight k and t be a non-constant automorphic function on X. It is known since the nineteen 
century that the functions F,tF, . .. ,r k F, as functions of t, are solutions of a certain 
(k + l)-st linear ordinary differential equations with algebraic functions as coefficients. 
(See lfTTiri"4l .) In particular, since i'(r) is a meromorphic automorphic form of weight 2, 
the function i'(r) 1 / 2 as a function of t, satisfies a second-order linear ordinary differential 
equation. We call this differential equation the Schwarzian differential equation associated 
to t, which has the following properties. 

Proposition 1 ([15, Proposition 6]). Let X be a Shimura curve of genus zero with elliptic 
points n, . . . ,T r of order e±, . . . , e r , respectively. Let t(r) be a Hauptmodul of X and set 
a L = t(ri), i — 1, . . . , r. Then t '(t) 1 / 2 ', as a function oft, satisfies the differential equation 

^F + Q { t)F = 0, 



where 



Q(t) = T 



1 ^ 1-1/e? , _ Bj 



4 ^ it- a.,) 2 ^ t- a, 
for some constants Bj. Moreover, if dj ^ oc for all j, then the constants Bj satisfy 

jZB 3 =Y. + \{l 1/e 2 )) = J2 UjBi + k(l - l/e 2 )) = 0. 



j—i j=i >• ' j=i 

Also, if a r ~ oo, then Bj satisfy 

X>, = 0, 2 Lb, + i(l - l/e 2 )) = \{l - l/el). 

We remark that when the Shimura curve has genus and precisely 3 elliptic points, 
the relations among Bj are enough to determine the constants Bj. This reflects the fact 
in classical analysis that a second-order Fuchsian differential equation with exactly three 
singularities is completely determined by the local exponents. When the Shimura curve 
has more than 3 elliptic points, the relations are not enough to determine Bj. In literature, 
we refer to this kind of situations by saying that accessory parameters exist. In order to 
determine the accessory parameters, one usually tries to find an explicit cover of Shimura 
curves and use it to determine the Schwarzian differential equations associated to the two 
curves simultaneously. 

Now one of the key observations in lfl5l is that the analytic behavior of the Hauptmodul 
t'(r) is very easy to determine and from this, one can work out a basis for the space of 
automorphic forms of even weight k in terms of t'{r). 

Proposition 2 (|15, Theorem 4]). Assume that a Shimura curve X has genus zero with 
elliptic points t%, . . . , r r of order e±, . . . , e r , respectively. Let t(r) be a Hauptmodul of X 
and set = t(r i ), i = 1, . . . , r. For a positive even integer k > 4, let 



d k = dim S k {X) = 1 - k + ^ 

3 = 1 



*/ _ JL 
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Then a basis for the space of automorphic forms of weight k on X is 

t'{r) k 'H{ry n (*(r)-ai)- Lfc(1 " 1/ei)/2J , j = 0, . . . , d k - 1. 

The combination of two propositions shows that all automorphic forms of a given even 
weight k can be expressed in terms of the solutions of the Schwarzian differential equa- 
tion. In [ 15 1, the author developed a method for computing Hecke operators relative to the 
basis in Proposition^ The key ingredients are the Jacquet-Langlands correspondence and 
explicit covers of Shimura curves. We refer the reader to [15 | for details. 

3. Computing modular equations for Shimura curves 

In this section, we will present a method for computing modular equations for Shimura 
curves, under the working assumptions that the Schwarzian differential equations have 
been determined and, for a certain prime fixed po, the matrices for the Hecke operator T po 
with respect to the bases in Proposition|2]have already been computed for sufficiently many 
k according to the recipe in ITT31 . 

Let X be a Shimura curve of genus of the form Xq(1)/W for some subgroup W 
of the group of Atkin-Lehner involutions and T be the discrete subgroup of SL(2,R) 

corresponding to X. For a prime p not dividing D, let jj = (^ a c 3 b d 3 ^j G GL + (2,M), 

j = 0, . . . ,p, be representatives for the cosets defining the Hecke operator T p , that is, 

jb/2-iy (det 7j ) fc / 2 

Then the modular equation Q p (x, y) = of level p is the polynomial relation between a 
Hauptmodul t(r) for X and t(jjr) for arbitrary j. 

Pick any nonzero automorphic form F(t) on X with the smallest possible weight k. 
For convenience, we let 

(det 7j ) fc / 2 

(CjT + dj) h 

the summands in ([T). Then any symmetric sum of Fj(r)/ F(t), j = 0, . . . ,p, will be an 
automorphic function on X and hence can be expressed as a rational function of £(t). In 
other words, there is a rational function ^(z, x) such that 

(3) =n(*-§S 

3=0 V V ; 

That is, we have 



(2) Fj(r)= r:? , J , k Fhjr), 



F(t) 
for all j. 

Lemma 3. Let ^{z, x) be the rational function defined by d3j. Let ^'(z, y) be the rational 
function such that 

which is a polynomial in z. Let R(x,y) be the resultant of^(z,x) and ty'(z,y) with 
respect to the variable z. Then the modular polynomial $ p (x, y) appears as one of the 
irreducible factors over C in the numerator of the rational function R(x, y). 
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Proof. Let j'j = (det 7j)t,- . There is an integer j', < j' < p, such that 7^ G r7j' . We 
apply the action of 7^ to (0]) and get 

*(w ,(7) ' T, )= a 



*' (' : §#.<(T-^))=0 



In other words, we have 

H/' ( 

F(r) 

From this equality and |@), we see that if we let R(x,y) be the resultant of ^5>(z,x) and 
^'(z, y) with respect to the variable z, then R(t(r), 4(7j't)) = 0. In particular, the mod- 
ular polynomial Q p (x, y) appears as an irreducible factor of R(x, y) over C. This proves 
the lemma. □ 

We remark that, in general, there will be more than one irreducible factors in the nu- 
merator of the resultant R(x, y), but it is not difficult to determine which one corresponds 
to $ p (x, y). For example, we can use the fact that the roots of x) are the coordi- 

nates of certain CM-points (also known as singular moduli) to test which factor of R(x, y) 
is Thus, from the above discussion, we see that the most critical part of the 

calculation is the determination of the rational function ^(z, x), which we address now. 

From Newton's identity, we know that the problem of determining the symmetric sums 
of Fj(r)/F(r) in (O is equivalent to that of determining (Fj/F) m for m = 1, . . . ,p + 1. 
Observe that the sum of F™ is equal to p 1 - fcm / 2 T p F, by the definition of the Hecke 
operator T p . Thus, to determine (Fj / F) m , one just needs to know how the Hecke operator 
T p acts on the basis given in Proposition|2] This is where the work of [ 15 1 comes into play. 

In 11151 , we developed a method to compute T po , p a prime not dividing D, for arbi- 
trary weight k with respect to the basis given in Proposition [2] assuming an explicit cover 
Xq(po)/W —> X®(l)/W is known. Now recall the following explicit version of the 
Jacquet-Langlands correspondence. 

Proposition 4 ([6, 8 1). Let D be discriminant of an indefinite quaternion algebra over Q. 
Let N be a positive integer relatively prime to D. For an Eichler order O = 0{D, N) of 
level {D, N) and a positive even integer, let Sk(Y{0)) denote the space of automorphic 
forms on the Shimura curve X^(N). Then 

S k (T(0)) ~ S?- a ™(DN) := Sr w (r a (dD))^ 

d\N m\N/d 

as Hecke modules. Here 

5r w (r (^)) [m] = {/(mr) : /(r) G {T o{dD))} 

and S^ cw (To(dD)) denotes the newform subspace of cusp forms of weight k on To(dD). 
In other words, for each Hecke eigenform /(r) in S®' new (T (DN)), there corresponds a 
Hecke eigenform /(r) in Sk(X(0)) that shares the same Hecke eigenvalues. Moreover, for 
a prime divisor p of D, if the Atkin-Lehner involution W p acts on f by W p f = e p f, then 

Wj=-e p f. 

For the situation under our consideration, we have 

(5) sr w (r OD),wo^(r), 

where the left-hand side denotes the subspace of the newform subspace of cusp forms of 
weight k on ro(D) that has eigenvalue —1 for all w q G W, q primes, and the right-hand 
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side denotes the space of automorphic forms of weight k on X. Thus, assuming the Hecke 
operator T Po on the space S^ cw (Tq(D), W) has no repeated eigenvalues, one can obtain 
the matrices for T p with respect to the bases in Proposition [2] from those for T po and the 
Fourier coefficients of Hecke eigenforms in S'^ cw (Tq(D), W). 

In summary, to compute the modular equation & p (x, y), we follow the following steps, 
assuming that the Schwarzian differential equation associated to X — Xq{1)/W and an 
explicit cover Xq (po)/W — > X^(l)/W are known for some prime po not dividing D. 

(a) Pick a nonzero automorphic form F on X of the smallest possible weight k, ex- 
pressed in the form given in Proposition!!] 

(b) Compute the matrices for T Pa with respect to the basis in Proposition|2]for weights 
k, 2k, . . . , (p + l)k using the method in |fT31 . 

(c) Compute Fourier coefficients of Hecke eigenforms in the space S^ cw (Tq(D), W) 
in ^ to the precision of p terms (using MAGMA ||2) or SAGE IfTUl ). 

(d) Compute the matrices for T p with respect to the basis in Proposition|2]using infor- 
mations from Steps (b) and (c). This gives us the expressions of Y^=o(Fj I 'F) m 
in terms of t(r), where Fj are defined by @. 

(e) Use Newton's identity to convert expressions for X^=o (^j I F) m t0 those for sym- 
metric sums of Fj/F and hence determine the rational function &(z, x) in (O. 

(f) Set ty'(z.y) = z p+1 ^>(l/ z,y). Compute and factorize the resultant R(x,y) of 
\l/(z, x) and ^'(z, y) with respect to the variable z. 

(h) Determine which irreducible factor of the numerator of R(x, y) is the modular 
polynomial ^ p (x, y) by using the fact that the roots of ^ p (x, x) are coordinates 
of some CM-points (also known as singular moduli) on X. 
We now work out an example in details. 

4. An example 

In this section, we will work out the modular equation of level 7 for the Shimura curve 
X = Xq°(1)/W w . (Note that this case was not covered in gj.) 

The curve X has 4 elliptic points of orders 2,2,2,3, which we denote by P2, P^, P2, P3, 
respectively. These are CM-points of discriminants —8, —20, —40, and —3, respectively. 
According to [4|, there is a Hauptmodul <(t) that takes values 00, 2, and 27, and 0, respec- 
tively. Using the covering X%° (3) /W w -> X 10 (l)/Wi , we |fT31 Equation (9)] found that 
the Schwarzian differential equation associated to t(r) is 

d 2 „ St 4 - lim 3 + 3157i 2 - 7296i + 10368 „ „ 

F H F = 0. 

dt 2 m 2 {t-2) 2 {t-27) 2 

Thus, by Proposition Q] near the point P3, the i-expansion of t'(r) is the square of a linear 
combination of two solutions 

„ , . , /, / 10 18539 9 168605 , 107269219465 , 

U; ~ V 81 839808 25509168 46548313473024 

„ , , 2/ o A 5 99095 2 8353325 , 851170821485 . 

Fo(t) = t ' 1 t t t t 

w V 81 5878656 1428513408 385081502367744 

of the differential equation above. Thus, by Proposition |2] a basis for the space of auto 

morphic forms of even weight k on X is 

(6) 

, t^{F 1 {t)-CF 2 {t)) k ' 2 
J*J < l*/3J(i _ t/2)Lfc/4J(l -</27)Lfe/4J ' 3 ^■■■' CLk 



k 




k 


_4_ 


+ 


3 
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for some complex number C. In particular, the one-dimensional space of automorphic 
forms of weight 4 on X is spanned by 



F(r) = h,Mr))- 



Let Fj(t), j — 0, ... ,7, be defined by (ffj with p = 7. According to the recipe described 
in Section [3] we first need to compute the matrices for T3 with respect to the basis in 
for weights 4m, m = 1, . . . , 8. This has already been done in lfT31 Appendix D]. We found 
that the matrices of the Hecke operators T3 are 



12 
16 
20 

24 
28 
32 



28 

468 ~98\ 
-1728 136 J 
1728 490 \ 
34560 -3572 / 

-2268 -2450\ 
-328320 35992 / 
/ 227772 -272244 14406 \ 

-388800 -258192 12250 
V 2985984 711936 -199556/ 
420552 949620 -72030 \ 
-933120 4479732 -61250 
104509440 31147200 -196568/ 
/ 29821932 -5456052 360150 
95084928 -48253536 306250 
\1803534336 -618444288 19290988/ 



That is, for the integer k in the table, we have 



fe,i 



fe,i 



= A k 

\fk.d k J \fk,d k 
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The next informations we need are the Fourier expansions of Hecke eigenforms in the 
space S£ cw (r (10), -1, -1) for k = 4, 8, . . . , 32. By MAGMA 0, they are 

q + 2q 2 -8q 3 + 4q 7 + ■ ■ ■ , 

q + 2 3 q 2 + 28q 3 + ■■■ + 104q 7 + • • • , 

q + 2 5 q 2 + aq 3 + ■■■ + (-177a + 60500)g 7 + ■ ■ ■ , 

q + 2 V + aq 3 + ■■■ + (423a - 102460)g 7 + • • • , 

q + 2 9 g 2 + aq 3 + ■■■ + (-417a + 48562100)q 7 + • • • , 

q + 2 1 Y + aq 3 + ■ ■ ■ + 7^(-a 2 + 156628a + 145233723936)g 7 + • • • , 

q + 2 13 q 2 + aq 3 + ■■■ + i(lla 2 - 17897672a - 75168751256976)g 7 + • • ■ , 
84 

q + 2 15 q 2 + aq 3 + ■ ■ ■ + — (a 2 - 14315428a + 407319502919904)g 7 + • • ■ , 
216 

respectively. Here each a is a root of the characteristic polynomial of T$ for the correspond- 
ing weight, and is different at each occurrence. From these, we deduce that the matrices 
for the Hecke operator TV with respect to the bases in (|6]l are 



12 
16 
20 

24 
28 
32 



Bk 

-4 
104 

-22336 17346\ 
305856 36428/ 
628484 207270 \ 
14618880 -1613416/ 
49507856 1021650 \ 
136909440 33553436/ 
/ -826476664 -2549118572 216037178\ 

-4554273600 -3184965196 546964950 
V 27509870592 52096359168 934073672/ 
/ -91564144564 113245670860 5617780910 \ 

438283664640 412736094176 -12502504350 
V15396166471680 -2162558865600 -111965170324/ 
/ 4631981436536 -203996300396 50284263050 
-11858411062656 12584751782372 97180354950 
V 18304356630528 78355740427776 4460404162424; 



Noticing that F m = fAm,d im f° r an Y positive integer m, we read from the matrices above 
that 

E | = 4 t % = ™ E^ = ^(305856t- + 36428), 

J=0 j=0 3=0 
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7 F 4 1 7 F 5 1 

^ = — (14618880i _1 - 1613416), ^ = ^(136909440i _1 + 33553436), 
7 F. 6 1 

J2~jk = - TT (27509870592t- 2 + 52096359168^ + 934073672), 

i=o 

7 J* 17 1 

^2p7 = 713 (15396166471680*" 2 - 2162558865600^ 1 - 111965170324), 



7 F. 8 1 

X/ = yi^(18304356630528i- 2 + 78355740427776^ 1 + 4460404162424). 



Then the rational function ^(z, x) in ([3]) is equal to 

*0, a;) = -^-2 (678223072849a; 2 z 8 + 387556041628a; 2 z 7 + 7909306972a; 2 z 6 
- (527663765132a; 2 + 4114130940864x)z 5 
+ (46199115214a; 2 - 5360751039168x)z 4 
+ (72916497220a; 2 - 2228082208128a;)z 3 
+ (90698975500a; 2 - 75419213184a; + 10905601867776)z 2 
+ (-72866748500a; 2 + 2516798571840a; + 9093300682752)z 

+ (13624725625a; 2 - 487484222400a; + 10905601867776)). 

We then set ^>'(z,y) — z 8 5'(l/ z,y) and compute the resultant R(x,y) of "$>(z,x) and 
^'{z, y) with respect to the variable z. The numerator of R(x, y) has two irreducible fac- 
tors. To determine which one corresponds to the modular equation $7(3;, y) of level 7, we 
use the fact that the roots of $7(3;, a;) should be the coordinates of CM-points of discrimi- 
nants -3, -20, -27, -35, -40, -52, -115, -180, and -280, and these coordinates are 
all rational numbers. In fact, the coordinates of these CM-points were given in Table 3 of 
[4|. Those obtained numerically in [4| were later verified by Errthum [5 1 using Borcherds 
forms. (In general, if & p (x, y) is the modular equation of level p for X, then the zeros 
of the polynomial & p (x, x) should be coordinates of CM-points of discriminants of the 
form (s 2 - 4p)// 2 , (4s 2 - 8p)// 2 , (25s 2 - 20p)// 2 , and (100s 2 - 40p)// 2 , subject to 
the condition that optimal embeddings of imaginary quadratic order of given discriminant 
into the maximal order in the quaternion algebra over Q of discriminant 10 exist.) That 
is, $7(3;, a;) should factor into a product of linear factors over Q. Indeed, exactly one of 
the two irreducible factors has this property. This determines &j(x,x). The equation of 
<&7(a;, a;) is given in Section|6] 

5. Modular equations for X$(l)/W 6 

In this section, we consider the Shimura curve X® (1)/Wg. The Hauptmodul t is chosen 
such that it takes values 0, 1, and 00 at the CM-points of discriminant —24, —4, and —3, 
respectively. For a prime p ^ 2, 3, we let & p (x, y) denote the modular equation of level 
p for the Shimura curve A^q(1)/W6- We have computed $ p (x, y) for primes up to 19, 
but because the coefficients are very big, here we only list the equation for p = 7. (The 
equation for p = 5, with a slight change of variables, is contained in Appendix A of lfT31l .) 
Files containing equations of other levels are available upon request. 
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We note that for p — 7, an explicit cover X§(7)/We — > X$(l)/We have aheady 
been determined by Elkies 101 . It is easier to use this explicit cover to obtain the modular 
equation. 

Write $7(2;, y) as a$(x)y 8 + • • • + a (x). Then 

a {x) = (262254607552729a; 2 - 121636570723920a; + 501956755356672) 2 , 



ai (a) = -2094804319407294388015256752187970421389152944241 1744000000a; 8 

- 972224700022233770983527054177244292949817108052713267466888a; 7 
+ 11273421679418251606098370957430742480477581392613807611381824a; 6 
+ 1668909333131109923039277853641492658428121008149001087888384a; 5 

- 63984186761131298882488446421883340458529413111323886398291968a 4 
+ 44207311496458731717072123962700099989994411322119128098799616a 3 
+ 16376343297097580660696811581414332932110848991772493800800256a; 2 

- 10266396414223579734319217958394551183170348131151337958146048a; 

- 61535122440667017884769009139405588503624662039298456944640, 



a 2 (a) = 38652996866188686573666243770903500386081949440000000000000a 8 

- 754695394662453117523773206823864267539586547407817536000000a; 7 

+ 106920403569574588553159546946652660456141317611471702546206108a; 6 

- 304712233218468979235116193697598666589103401841722665615811904a; 5 
+ 760829927783454753461831780276599490081935198148332681198238208a; 4 

- 838656978040090256476182485622824884772400584504768719373500416a 3 
+ 279281854044007197025412087758954247401641295020911801816973312a; 2 
+ 16376343297097580660696811581414332932110848991772493800800256a 
+ 155040070253203796103748431879063330519757479108093010771968, 



a s (a) = -38038344753469448835190136766681303367680000000000000000000a 8 
+ 5356342063300144660882808541232703156528549197824000000000000a 7 
+ 76321377218916949783840090111683558976051518322915065600000000a 6 

- 474977723264104825683894637831760410975632622470023372974366136a 5 
+ 167509988385063268988566748828192523079606899212802760635355200a; 4 
+ 929561753522872554487774286564540945788659043162410938085107712a 3 

- 838656978040090256476182485622824884772400584504768719373500416a 2 
+ 44207311496458731717072123962700099989994411322119128098799616a 

- 100063182095188848153455102416032183737981837866201188925440, 
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a 4 (x) = 21056325026580996459284526662615040000000000000000000000000a; 8 

- 2676503627950279289618572142817040185729024000000000000000000a; 7 
+ 74717592843661815288079677324740635068097597899264000000000000a; 6 
+ 99843100080788043313984402843762031194115908069343957280000000a; 5 

- 643716270399690372725464450317965610492381617103712363036843130a; 4 
+ 167509988385063268988566748828192523079606899212802760635355200a; 3 
+ 760829927783454753461831780276599490081935198148332681198238208a; 2 

- 63984186761131298882488446421883340458529413111323886398291968a; 
+ 127566357194950739797651934018366979046647830016645952241664, 

a 5 (as) = -6216447485624191385037963264000000000000000000000000000000a; 8 

- 969897592146461089985813356429049856000000000000000000000000a; 7 

- 4716808719470806728440581205159667124740096000000000000000000a; 6 

+ 100946931718713025406590060876901549085887873170944000000000000a; 5 
+ 99843100080788043313984402843762031194115908069343957280000000a; 4 
-474977723264104825683894637831760410975632622470023372974366136a; 3 

- 304712233218468979235116193697598666589103401841722665615811904a; 2 
+ 1668909333131109923039277853641492658428121008149001087888384a; 

- 52279464855900563088298143909713657632323089372536352686080, 
a 6 (x) = 764699349893278334976000000000000000000000000000000000000a; 8 

- 15347209366323527757160513536000000000000000000000000000000a; 7 

- 1253143615853758839349272867028205568000000000000000000000000a; 6 

- 4716808719470806728440581205159667124740096000000000000000000a; 5 
+ 74717592843661815288079677324740635068097597899264000000000000a; 4 
+ 76321377218916949783840090111683558976051518322915065600000000a; 3 
+ 106920403569574588553159546946652660456141317611471702546206108a; 2 
+ 11273421679418251606098370957430742480477581392613807611381824a; 
+ 42321163987340612844969199376648358822165115192172887352832, 

a 7 (x) = 997738013984972341248000000000000000000000000000000000000a; 7 

- 15347209366323527757160513536000000000000000000000000000000a; 6 

- 969897592146461089985813356429049856000000000000000000000000a; 5 

- 2676503627950279289618572142817040185729024000000000000000000a; 4 
+ 5356342063300144660882808541232703156528549197824000000000000a; 3 

- 754695394662453117523773206823864267539586547407817536000000a; 2 

- 972224700022233770983527054177244292949817108052713267466888a; 

- 8775937874145070526339174290824946318421255782225802867520, 

a s (x) = (3024000000x - 4097152081) 6 . 
We next give a short list of irrational singular moduli obtained by factorizing Q p (x, x) 
forp < 19. Note that the norms of these coordinates were already computed in (5). Our 
computation yields their exact values, not just their norms. (Note that all the rational 
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singular moduli for Xq(1)/Wq were numerically determined in [4| and later verified in 
0.) 



coordinates 



-264 
-276 
-136 
-195 
-420 
-456 
-219 



26198073a: 2 - 53148528a; + 42719296 
1771561a: 2 + 328736070x - 206586207 
8703679193a: 2 + 9087757328a; + 52236012608 
1048576000000a: 2 - 1812440448000a; + 884218628241 
377149515625a; 2 - 1357411709250a; + 3732357388761 
42761175875209a: 2 - 43342656704496a; + 31050765799488 
1320903770112a; 2 + 125096383466496a; - 128030514913249 



2 f, 19 2 43 2 /3 a ll J 

3 8 23 1 37 2 /11 6 

2 6 13 4 41 2 /11 6 17 2 

3 s 13 2 19 2 47 2 /2 26 5 6 

3 8 23 2 61 2 /5 6 17 4 

2 6 3 9 19 1 67 2 /11 6 17 4 

13 4 23 2 41 2 71 2 /2 26 3 9 



6. Modular equations for Xq°(1)/W w 

In this section, we consider the Shimura curve Xq°(1)/Wiq. The Hauptmodul t is 
chosen such that it takes values 0, 2, 27, and oo at the CM-points of discriminant —3, —20, 
—40, and —8, respectively. For a prime p ^ 2,5, we let y) denote the modular 

equation of level p for the Shimura curve Xq°(1) /Wiq. We have computed & p (x, y) for 
primes up to 23, but here we only list the equations for p = 7. Files containing equations 
of other levels are available upon request. (The equation for p = 3 is given in Section 5 of 

ma.) 

Write $7(3;, y) as as(x)y 8 + • • • + ao(x). Then 

a {x) = a; 2 (278055625x 2 - 9948657600a; + 222563303424) 3 , 

ai(x) = -22595841233117014160156250a; 8 
+ 28984697923550260459599062500a; 7 

- 1645302326393970545433073200000a; 6 
+ 6642716791333154856889291776000a; 5 
+ 602202478597261 1 11037649172561920a; 4 

- 4910763943297233363215784691630080a; 3 

- 8496599535482361705366127525232640a; 2 
+ 9192477793078594137923003873230848a;, 

o 2 (a:) = 227107188088738387105484375a; 8 

- 10394755051722304542380018750a; 7 
+ 2778503946540202482367606147300a; 6 

- 87342791194115920501573603009536a; 5 
+ 509955665442187448278395883376640a; 4 
+ 4112275508763979970115841213071360a; 3 

- 2548787139523595922308471728373760a; 2 

- 8496599535482361705366127525232640a; 
+ 11024545045545306955620945204609024, 
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a 3 (x) = -150467239701697875855937500a; 8 

- 66608231302181029579976003090a; 7 
+ 1779457212256262244783949440310a- 6 
+ 2537754854567828839953906813980a; 5 

- 567333242539975544914354048540992a; 4 

- 239926105649803444132294072811520a; 3 
+ 4112275508763979970115841213071360a; 2 

- 4910763943297233363215784691630080a; 

- 1478403072292996522939287129292800, 

a 4 {x) = 591743143857926286378268151a; 8 
+ 3202866389133529516632265220x 7 

- 35406367266732536603 1 194372695a; 6 
+ 22584679924916438377343486492770a; 5 
+ 91353706850814722447426682881030x 4 

- 567333242539975544914354048540992a; 3 
+ 509955665442187448278395883376640a; 2 
+ 602202478597261111037649172561920x 
+ 107405110735755767009079459840000, 

o 5 (a;) = -178678323923978351074342650a; 8 
+ 71993226806139685483910592a- 7 

- 382859191369766872635900582760a; 6 

- 3892236690880311638995825770460a; 5 
+ 22584679924916438377343486492770x 4 
+ 2537754854567828839953906813980a; 3 

- 87342791194115920501573603009536a; 2 
+ 6642716791333154856889291776000a- 

- 4678710992156906294833152000000, 

a 6 (x) = 38853800265463031036085625a; 8 
+ 2535839489102333226932059850a- 7 
+ 55519696363179781200215763218a; 6 

- 382859191369766872635900582760a; 5 

- 354063672667325366031 194372695a; 4 
+ 1779457212256262244783949440310a; 3 
+ 2778503946540202482367606147300a; 2 

- 1645302326393970545433073200000a; 
+ 134184722882776668446400000000, 
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a 7 (x) = -3904463975228746224750000a: 8 

- 279343739152597734767051250a; 7 
+ 2535839489102333226932059850a: 6 
+ 71993226806139685483910592a; 5 
+ 3202866389133529516632265220a: 4 

- 66608231302181029579976003090a; 3 

- 10394755051722304542380018750a; 2 
+ 28984697923550260459599062500a; 

- 2307539315546608933125000000, 

a 8 {x) = (528806915000a; 4 - 3691767131325a; 3 + 23850535042201a; 2 

- 2436693984375a; + 4636577546875) 2 . 

Here we give a short table of irrational singular moduli obtained by factorizing Q p (x, x) 
forp < 23. Note that the norms of these singular moduli were already computed in [5] us- 
ing Borcherds forms. (All the rational singular moduli for Xq ( 1 ) / Wiq were numerically 
determined in [4| and later verified in ||5l.) 



d 


coordinates 


norms 


-68 


25a; 2 - 164a; + 500 


2 2 5 x 


-260 


31213a; 2 - 7252a; + 19652 


2 2 17 3 /7 4 13 1 


-360 


3844a; 2 + 16299a; + 132651 


3 3 17 3 /2 2 31 2 


-152 


1250a; 3 - 62475a; 2 + 64048a; - 166375 


11 3 /2 1 5 1 


-132 


4225a; 2 - 27108a; + 364500 


2 2 3 6 5 1 /13 2 


-168 


122500a; 2 - 647325a; + 970299 


3 6 ll 3 /2 2 5 4 7 2 


-195 


142129a; 2 + 12489984a; - 2985984 


2 12 3 6 /13 2 29 2 


-155 


74431a; 2 - 1204224a; - 5451776 


2 12 n 3 /7 4 31 l 


-440 


114244a; 3 - 13732264a; 2 + 20581253a; - 16194277 


11 3 23 3 /2 2 13 4 


-83 


105625a; 3 - 2894400a; 2 + 30224384a; + 32768000 


2 18 /5 1 13 2 


-420 


67081a; 2 + 431676a; + 71118324 


2 2 3 6 29 3 /7 2 37 2 


-580 


13766909a; 2 - 22395636a; + 200973636 


2 2 3 6 41 3 /13 2 29 1 53 2 


-315 


72812089a; 2 + 88682944a; + 398688256 


2 15 23 1 /7 2 53 2 


-820 


3158777209a; 2 - 19708111836a; - 598885164 


2 2 3 6 59 3 /7 4 13 4 37 2 


-355 


253887551a; 2 - 6128728704a; + 3974344704 


2 12 3 6 11 3 /31 2 61 2 71 


-660 


4726150009a; 2 - 10741876236a; + 8174192436 


2 2 3 9 47 3 /7 4 23 2 61 2 


-680 


995781136a; 3 - 8234027816a; 2 + 166339878257a; - 450688409963 


ll 3 17 3 41 3 /2 4 7 6 23 2 


-435 
-920 


4460503369a; 2 - 298408090176a; + 25350000869376 
287478606136k 5 -3389457456860a: 4 -4652502434310x 3 
+ 2 16562501495k 2 + 50275785584380a: - 60944362239977 


2 18 3 9 17 3 /7 4 29 2 47 2 
23 2 59 3 /2 3 29 1 47 2 
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